Strongly enhanced light-matter interaction in a hybrid photonic-plasmonic resonator 
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We propose a hybrid photonic-plasmonic resonant structure which consists of a metal nanoparticle 
(MNP) and a whispering gallery mode (WGM) microcavity. It is found that the hybrid mode enables 
a strong interaction between the light and matter, and the single-atom cooperativity is enhanced by 
more than two orders of magnitude compared to that in a bare WGM microcavity. This remarkable 
improvement originates from two aspects: (1) the MNP offers a highly enhanced local field in the 
vicinity of an emitter, and (2), surprisingly, the high-Q property of WGMs can be maintained in 
the presence of the MNP. Thus the present system has great advantages over a single microcavity 
or a single MNP, and holds great potential in quantum optics, nonlinear optics and highly sensitive 
biosening. 

PACS numbers: 42.50.Pq, 42.50.Ct, 42.50.Dv, 78.67.-n 



Owing to the size mismatch between light and single 
emitters such as single atoms, the interaction between 
them is very weak, so that it is of importance to create a 
light-matter interface enabling strong interactions. One 
way to bridge this mismatch is to employ the strong in- 
teraction within cavity quantum electrodynamics (QED) 
PjQ. Cavity QED offers an almost ideal platform for the 
study of physics at the interface of classical and quantum 
mechanics, and provides a technology for various devices 
in the field of quantum information 3-5]. Experiments 
on strong coupling regime in cavity QED have made great 
advances over the past two decades @. Among them, 
whispering gallery mode (WGM) microcavities 0) are 
promising because they possess ultrahigh quality (Q) fac- 
tor and allow for mass production on a chip. However, 
the relatively large cavity mode volume makes it difficult 
to realize strong coupling. On the other hand, due to the 
localized surface plasmon resonance (LSPR) Q, metal 
nanoparticles (MNPs) [§] enable subwavelength confine- 
ment of the optical field [Tol - [l3| . Unfortunately, MNPs 
suffer from serious absorption and scattering losses. 

Against this backdrop, in this Letter, taking advan- 
tages from both ultralow-loss WGMs and highly localized 
plasmon, we propose a WGM microcavity-MNP resonant 
system. In this composite system, the high-Q WGM 
microcavity serves as a low-loss storage of the optical 
field, while the MNP plays the role of an optical antenna 
which creates a hot spot and magnifies the local optical 
field. Remarkably, the high-Q property of WGMs can be 
maintained in the presence of the MNP. As a result, the 
cooperativity parameter (defined as C — 2G 2 /kj s [lij ]. 
with G being the single photon coupling strength, j s the 
spontaneous decay of the emitter and n the decay of the 
cavity field) achieves a more than 100-fold increase com- 
pared with that of the WGM cavity alone. It should be 
pointed out that, this composite cavity QED structure is 
significantly different from previous designs where a silica 
disk or toroid was completely covered with a metal lay er, 
which led to strong degrading of the Q- factor lfj, [l7| • 




FIG. 1: (Color online) (a) Sketch of the MNP- WGM compos- 
ite system (not to scale), (b) Zoom- in view of the MNP and 
the dipole emitter. The inset is the energy diagram of the 
dipole. (c)-(e) Three-dimensional finite element simulations 
of the electric field distribution. Top view of the electric field 
profile without (c), and with (d) the MNP. (e) Zoom-in view of 
the field distribution near the MNP (azimuthal cross-section). 
The red arrows show the electric field directions. The vertical 
line on the left of the MNP is the cavity boundary. The col- 
ors in (c)-(e) represent the norm of the electric field (|.E7|). (f) 
Comparison of the electric field with (blue dot) and without 
(light gray dot) MNP for 9 = (left) and 6 = vr/2 (right). 
The red solid curves are the theoretical result (see text). To 
save simulation resources, here we use a relatively small mi- 
crocavity (major radius 0.8 /im) and a relatively large MNP 
(r m = 30 nm). 



Figure Q] illustrates a schematic of the system. A 
MNP is located onto the surface of a microtoroidal cav- 
ity [IH which supports twin counter-propagating WGMs 
with degenerate frequency lu c (neglecting the coupling 
between these two modes arising from surface rough- 
nesses induced scattering). A dipole emitter (ground 
state \g), excited state |e), energy spacing huj e and dipole 
moment fj,e x ) is placed in the vicinity of the MNP. Ex- 
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perimentally we can use atomic force microscope (AFM) 
manipulation to controllably position the MNP and the 
dipole emitter [2(j. Alternatively, we can also use 
gold-coated fiber taper tip to play the role of the MNP, 
in analogy to the method utilizing fiber tip to simulate 
dielectric nanosphere [2l|. A tapered fiber is used to 
couple light into and out of the microcavity. We con- 
sider TE-polarized WGMs, with the dominant electric 
field component in the e x direction. In the absence 
of the MNP, the Jaynes-Cummings Hamiltonian reads 



H, 



4= 



\uj c a z +G c I] n (a T t l 0-_ +a„er + ), where 
(J-=(J+= \g) (e\, a z = |e) (e| - \g) (g\, a n (a*) denotes an- 
nihilation (creation) operator of the WGMs with n — 
CW (CCW) corresponding to the clockwise (counter- 
clockwise) propagating mode. The single photon cou- 
pling strength (vacuum Rabi frequency) is given by G c = 

1 /2 

^/ c (R)(cj c / (2he £ c V c )) 1 , where £o is the permittivity of 
vacuum, e c is the relative permittivity of the microcavity, 
V c and / C (R) denote the mode volume and normalized 
field distribution of the WGMs, respectively. 

Now we focus on the case that a spherical MNP is 
placed in the vicinity of the microcavity and the dipole 
emitter, with the position R = Rq. The radius of the 
MNP r m is much smaller than the light wavelength so 
that the interaction arising from the MNP is governed 
by electrostatics rather than electrodynamics [22[ , which 
is confirmed by our numerical simulations (Fig. I (c)-(f)). 
By solving Laplace's equation, we obtain the positive fre- 
quency component of the total field inside and outside 



the MNP: 



E<+2(r,0) 



(1 - (3)E { c +, (Tl )e x ,(r < r m ), 
E ( c +) (R )e x +/3(rfJr 3 )E { c +) (R ) (2e r cos 6+ 
egsin9),(r > r m ). Here the dipole emitter-induced 
field has been omitted in the weak excitation limit; 
r = |R — Ro| and 9 form a polar coordinates system; the 
complex coefficient (3 is given by (3 — (e m — £b)/(e m +2£b), 
where e m (eh) denotes the relative permittivity of the 
MNP (the surrounding medium) . We can use the Drude 
dispersion relation £ m (co) = I — bj^/[uj{bj + i~f m )}, where 
lu p is the bulk plasma frequency and 7, n is the damping 
rate which accounts for energy dissipation due to ohmic 
losses. The LSPR [|| occurs with the dipolar plasmonic 
resonance frequency w sp satisfying Re[e m (u; S p)] = — 2eb. 
In this case, |/3| ^ I, and the local field in the vicinity of 
the MNP is enhanced significantly. 

The response of the MNP results in the re- 
distribution of the cavity mode field, forming a hy- 
brid photonic-plasmonic mode. The modified maximum 
field strength reads |£? c ,m,max| = /c(Ro) 1(1 + 2^)E c , max |, 
and thereby the mode volume decreases to be V^ m = 
£ c V c /[et, \l + 2/3 1 2 / 2 (Ro)], which approximately scales as 
times smaller than bare cavity case. Thus the MNP 
effectively plays the role of an optical antenna which (i) 
confines the photon energy, (ii) reduces the mode vol- 
ume, and (iii) magnifies the local field. For the longitu- 
dinal case (6 = 0), the modified single photon coupling 



strength G Cim is calculated to be |l + 2/3r 3 n /r 3 | times as 
large as G c . 

In this composite cavity QED system, a key question 
is whether the MNP significantly degrades the high-Q 
property of WGMs. The MNP-induced decay includes 
two contributions: scattering and absorption losses. For 
the subwavelength MNP, the scattering interaction can 
be modeled in a dipole approximation, where the elec- 
tric field of the input wave induces a dipole moment 
in it 21]. On one hand, this scattering interaction 
mixes the two counter-propagating modes, described by 
the Hamiltonian H\ = — |p m • E c m (R ), where p m = 

£Q£hCtEc + \Ho)e x + H.c. is the polarization and a = 
47rr^ 1 /3 is the polarizability of the spherical MNP. It can 
be simplified as Hi = hJ2 n „' a n a n', where n, n' = CW, 
CCW, and h = 27rrf n e b c 1 ; c \(3\ 2 / c 2 (R )/(e c V c ) is the cou- 
pling strength. On the other hand, the scattering results 
in the decay from WGMs to reservoir modes, described 
by the Hamiltonian H 2 = — ^Y^jiPm ' E?(Ro) + Pj • 
E C)1D (R )), where Ej(R ) is the electric field for the jih 
reservoir mode and Pj is reservoir-induced polarization of 
the MNP. Taking the Weisskopf-Wigner semi- QED treat- 
ment [2lj, the scattering process results in radiation en- 
ergy decay of the WGMs, which takes an analogous form 
of spontaneous emission. After tracing out the freedom 
of the reservoir modes, we obtain the effective Hamilto- 
nian i?2,eff = —ii^K J2 n n' a n a m with the damping rate 
KR = e b /2 (47rri) 2 |/3| 4 Wc 4 / c 2 (Ro)/(6^c 3 e c K), where c is 
the speed of light in vacuum. 

The absorption of the MNP results in ohmic losses, 
characterized by the imaginary part of the metal's 
permittivity. Thus this energy dissipation can be 
described by the non-Hermitian Hamiltonian i?3 iC ff = 
5 £ J|R-R |<r' m Im M( we m(^))/w] | £'c,m(R) 2 rfR 3 , 
re-written as H^^s = ~in m J2 n n' a n a nh with the decay 
rate n m = Anr^ |1 - /3| 2 w2 7m / c 2 (R )/(3e c w c 2 y c ). In 
current case, the energy fraction of the hybrid mode 
within the MNP is very small (< 0.01%), resulting a 
minor K m . 

Now we arrive at the overall Hamiltonian of the com- 
posite system H c m = H c + Hi + H 2 , c t[ + H 3 c g. Next 
we consider specific examples. Different from low-Q pho- 
tonic crystal nanocavity [23j, here we focus on the en- 
hanced light-matter interaction in a silica microtoroidal 
cavity [18| which possesses a relatively high intrinsic qual- 
ity factor Qo = I0 7 . We use a typical toroidal microcav- 
ity with the major and minor radii being 30 /xm, 3 /im, 
respectively. For such a cavity in air, e c = 1.45 2 , £b = 1, 
we obtain V c ~ 200 /im 3 (corresponding to 4000(A/n) 3 ) 
and / c (Ro) ~ 0.3 using finite element simulations. We 
set the cavity-taper coupling strength Ki = 5kq, where 
kq = uj c /Qo denotes the intrinsic damping of the WGMs. 
We use large Ki here so as to realize near critical cou- 
pling in the presence of the MNP which brings about 
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FIG. 2: (Color online) (a) Contour plot of the cooperativity 
enhancement (C c , m /C c ) vs. r m and d. (b) Parameters {G c , 
G c ,m, h, Ko, kr, K m }/27r vs. r m . (c) Cooperativity with the 
presence of the MNP (C c . m ) and with its absence (C c ). The 
horizontal line in (a) corresponds to d = 3 nm. The vertical 
lines in (a)-(c) corresponds to r m = 12 nm. (d) Transmission 
spectra for r m = 12 nm and d = 3 nm in the presence of the 
dipole, with (solid line) and without the MNP (dashed line). 
{G c , G c ,m, h, k , kr, K m }/27r ={760, 9000, 170, 55, 80, 30} 
MHz. The input-cavity detuning A = uj — ui c . The inset is 
for no dipole case. For (b)-(d), d = 3 nm; For (a)-(d), the 
cavity-MNP detuning A sp = 0. 



additional decays. For a gold MNP, the permittivity 
can be extracted from the experimental data .24], with 
w p ~ 6 x 10 15 Hz (corresponding to the LSPR wave- 
length of 540 nm) and 7, n ~ 3 x 10 14 Hz. For the 
dipole emitter, we can use chemically synthesized cad- 
mium selenide (CdSe) quantum dot, with the dipole mo- 
ment n = 2.4 x 10~ 28 C-m [l3j], the spontaneous emission 
rate 7 S = ey 2 / u 2 w|/(37re Sc 3 ) = 27r x 1.6 GHz and the 
emitter-cavity detuning A ec = lu c — lu c = 0. To focus on 
the physics, the size of the quantum dot is assumed to 
be small so that its mesoscopic effects [25| can be ne- 
glected. Note that for a more general emitter like an 
atom or a dye molecule, the size can be considered nil. 
In Fig. [2fa) we plot the cooperativity enhancement as 
a function of r m and d (the distance between the dipole 
emitter and the surface of the MNP). It is shown that 
the enhancement can exceed more than two orders of 
magnitude, and there is enough parameter space for re- 
markable enhancement. As an example, we set d = 3 
nm, corresponding to the horizontal line in Fig. At 
this distance, the charge carrier tunneling can be s afely 
avoided [H,^. With a MNP of specific geometry [28j . 
the dipole's emission to dark multipolar plasmonic res- 
onances, which can be regarded as an addition sponta- 
neous decay of the dipole emitter and largely related to 
the concept of the "quenching" in the fluorescence mea- 
surement, will be significantly reduced to a minor effi- 
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FIG. 3: (Color online) (a) Cooperativity enhancement 
(Cc,m/C c ) vs. cavity-MNP detuning A sp /7 m for several typ- 
ical r m . (b) The cooperativity parameters C c , m , C\ >m (near- 
resonance, [Asp I <C 7m) and (7" m (off-resonance, |A sp | > j m ) 
for r m = 20 nm. 



ciency of 1% — 3%. The metal absorption in the hybrid 
mode, also related to "quenching" , is extremely weak 
due to the ultra-small energy fraction (< 0.01%) within 
the MNP, which differs from the pure plasmonic reso- 
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nant system. When LSPR occurs (A sp 
we obtain |/3| = 11.5. Figure E][b) plots the characteris- 
tic parameters depending on the radius of the MNP. It 
explicitly shows that the dipole-field single photon cou- 
pling strength with the help of the MNP (G Cjm , solid line) 
can be more than 10 times larger than that without the 
MNP (G c , dashed line), with achievable G Cim /27r ~ 13 
GHz. When r m is small, the extra decays (kr, « m ) are 
smaller than the intrinsic cavity decay kq- Therefore, for 
a small MNP, the enhancement of the local field can be 
quite strong, while the high-Q property of the system is 
surprisingly maintained. In Fig. [2jc), we plot the co- 
operativity parameter as a function of r m both in the 
presence of the MNP and in its absence. We can see that 
for r m ~ 12 nm, the cooperativity is enhanced more than 
100-folds, with the corresponding G c , m /27r = 9 GHz. 

To probe the enhanced light-matter interaction, we use 
weak monochromatic clockwise propagating light input 
aJ.'L. By utilizing the Heisenberg-Langevin equations and 
the input-output theory [2^ . [30 |. we obtain the cavity 
transmission spectra for r m = 12 nm in Fig. [UJd). When 
the MNP is absent, the vacuum Rabi splitting cannot be 
resolved. On the contrary, with the help of the MNP, a 
mode splitting of 2v2G c , m can be observed. This split- 
ting is completely the cavity QED effect and the scatter- 
ing of the MNP has no contribution to it. To demonstrate 
this, in the inset of Fig. [2jd), we plot the transmission 
spectra for the no dipole emitter case. As expected, no 
splitting is observed, where the asymmetry spectrum in 
the presence of the MNP is due to the scattering-induced 
mixing of the WGMs. 

One distinct property of LSPR is that the plasmon 
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FIG. 4: (Color online) Illustration of the fast quantum gate 
operation (a) and NOON state generation (b). 

excitation covers a broad bandwidth. In fact, the op- 
timal condition for the improvement of strong coupling 
does not occur at the exact cavity-MNP resonance for 
large-size MNP where the MNP-induced decay exceeds 
the original decay. This is verified in Fig. [3][a), which 
displays the cooperativity parameters for cavity-MNP de- 
tuning (A sp = lu c — w sp ) varying from —67m to 47 m . It is 
interesting that C c ,m/C c arrives at a local minimum at 
A sp = for a large r m . The underlying physics is that 
the MNP-induced decay contributes more than the field 
enhancement. To explain this more clearly, we plot the 
cooperativity for r m = 20 nm in Fig. [3Jb). In the near- 
resonance region, the loss induced by the MNP is domi- 
nant, and thereby C Cjm ~ C l cm = 2G^ m /[%(n R + K m )]; 
while in the off-resonance region, these decays decrease to 
a low level, yielding C c , m ~ C^ m = 2G 2 cm /[^ s (K + K%)]. 
It is shown that, for red cavity-MNP detuning case, the 
cooperativity can maintain a high value over a broad 
range. Note that A sp = 0, —67m corresponds to the 
wavelength 540 nm, 1130 nm, respectively. Moreover, 
even for r m = 30 nm, the cooperativity can still be en- 
hanced by 30-folds for a suitable detuning. 

Enhancing the strong coupling in cavity QED is of cru- 
cial importance in high-speed operation of quantum gate 
and fast generation of entangled state, as illustrated in 
Fig. QJa). These can be obtained among two or more 
identical MNP-atom molecules localized near the micro- 
cavity. The atoms interact with a shared cavity mode and 
a classical laser field is applied to control the evolution 
among them. Typically, to suppress the atomic spon- 
taneous emission, the system usually works in the large 
atom- field detuning limit (A oc ^> G) 0, [3lj] . For conven- 
tional cavity QED setup without MNPs, this large detun- 
ing limit is theoretically feasible, but may not be applica- 
ble because the effective interacting strength (oc G 2 /A ec ) 
decreases significantly and the Rabi oscillation becomes 
slow. On the contrary, with the help of MNPs, the effec- 
tive interaction is significantly enhanced. As an exam- 
ple, the effective interacting strength can be improved 
100-folds with the help of MNPs. This means that the 
quantum gate operation and the entanglement genera- 
tion can finish in a much shorter time, and the system 
allows for much more operations within the coherence 
time. 



Other examples include the fast generation of spin 
squeezed state [32j and N dipole emitters maximally 
entangled state (also known as NOON state) [33j]. To 
this end, N identical MNP-atom molecules are posi- 
tioned around the cavity. In the large detuning case, 
the effective Hamiltonian is given by H e g — x(Jz ~ J^), 
where x = 2G* m /(A ec - Nh) and J z = ±£f =1 <#. 
This Hamiltonian is the generalized one-axis-twisting 
type [34] ]. which implies that it is capable of perform- 
ing spin squeezing. To create the NOON state, dipole 
emitters are initially prepared in the superposition state 
2 ~ Ar/2 (l5) + |e))' 8Ar using a fast tt/2 pulse. Next the sys- 
tem undergoes free evolution for \t — Then a sec- 
ond 7r/2 pulse is applied to finally create the NOON state 
(\gf N + \ef N )/V2, as described in Fig. 1(b). Previ- 
ously this scheme is subjected to the weak nonlinearity 
and experimentally it is a challenge to increase the non- 
linearity so as to overcome the decoherence. In our pro- 
posed composite system the nonlinear coefficient x can 
be large enough so that the squeezing and the NOON 
state can be obtained in a short time, which is very im- 
portant in view of decoherence. Moreover, the additional 
dissipation induced by the MNP has little effect on the 
system because here the dipole emitters interact disper- 
sively with the cavity modes. 

In summary, we have studied a composite MNP-WGM 
coupling system. With the help of the MNP, the sin- 
gle photon coupling strength between a dipole emitter 
and the cavity system can be substantially enhanced, 
while the high-Q property is surprisingly preserved. The 
single-atom cooperativity obtains a more than two orders 
of magnitude increase compared with the conventional 
WGM cavity QED system. Also, the MNP-induced cou- 
pling enhancement covers a broad band up to several 
hundred nanometers so that the microcavity and the 
MNP do not have to be on resonance. The system is ap- 
plicable ranging from quantum optics to quantum infor- 
mation science. It also offers new opportunities to exploit 
both theoretical and experimental physics in stronger 
light-matter interaction regime, such as nonlinear op- 
tics (e.g., enhanced Raman and Rayleigh scattering) and 
highly sensitive biosening (e.g., single nanoparticle sizing 
with an improved detection limit). 
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